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Abstract 

Intrinsic aberrations are those which occur due to the finite length 
of the desired field configuration. They are often loosely ascribed to 
the fringing field. This is misleading as it implies that the effects can 
be minimized by shaping the fields. In fact, there is an irreducible 
component related to the broken symmetry. It is present even in the 
hard-edge limit, and moreover, the other (soft-edge) effects can be 
simply ascribed to the intrinsic aberration spread over a finite length. 

We rederive the aberration formulas for quadrupoles using a Hamil¬ 
tonian formalism. This allows for an easy comparison of electrostatic 
and magnetic quadrupoles. For different combinations of large and 
small emittances in the two transverse planes, it is found that in some 
situations electrostatic quadrupoles have lower aberrations, while in 
others, magnetic quadrupoles are better. As well, we discuss the ways 
in which existing transport codes handle quadrupole fringe fields. Pit- 
falls are pointed out and improvements proposed. 


1 Introduction 

A common prescription for quadrupole design is that the beam occupy no 
more than a certain fraction, usually 1/2 to 3/4, of the aperture. Another 
common prescription is that the longer the quadrupole is compared with its 
bore diameter, the better. As well, it is common practice to carefully round 
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the ends of the poles. We show that none of these common practices can be 
validated by the quadrupole dynamics np to 3 rd order in force. 

The misconceptions are perpetuated by existing transport codes like GI0S[1] 
and C0SY[2], because they allow one to make calculations with no fringe 
fields. Then when fringe fields are included, aberrations increase. In fact, 
the no-fringe-field cases are non-physical and such calculations should not be 
permitted in transport codes. The aberrations in question are not caused by 
the fringe fields, but by the broken symmetry inherent in a quadrupole of 
finite length. 

We start with the quadrupole Hamiltonian and find canonical transfor¬ 
mations for both the electrostatic and magnetic cases which eliminate the 
derivatives of the quadrupole strength up to 4 th order. In this way, we eas¬ 
ily reproduce the known aberration formulas, but with additional physical 
insight. 

2 Theory 

2.1 Electrostatic 

When using the longitudinal position as independent variable, the Hamilto¬ 
nian H is just the longitudinal momentum: 

H = -\Jp 2 o - 2 ~ P 2 X ~ P 2 y . (1) 

The electrostatic potential is <&(x,y,z) and po is the reference momentum. 
We will benefit from cleaner and more transparent notation if momenta are 
measured in units of p 0 . This has the additional benefit that to first order, 
p x = x'. Additionally, we let the potential <f> in units of the reference kinetic 
energy pg/(2m). Then 

l ( 2 ) 

We expand the square root to 4th order in coordinates and ignore the con¬ 
stant: 

^+Pl+P 2 y) +^($+pl+p 2 y) 2 . (3) 

To the same order, Laplace’s equation gives for the expansion of the 
quadrupole potential: 

^ = V(z)(x 2 -y 2 )-^-^(x 4 -y 4 ). (4) 
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The final Hamiltonian, correct to 4 th order is 
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The trouble with applying this to simple cases like thin lenses and hard- 
edge limits is the presence of V"(z), which becomes singular in those limits. 
In most cases, one sacrifices physical insight and simply traces particles with 
this Hamiltonian, using a more-or-less realistic function V(z). For example, 
the approach taken in GIOSjTj is to leave it up to the user to specify ‘fringe 
field integrals’ such as / V 2 dz through the fringe fields. However, this leaves 
much room for error; different integrals may not be realistic or consistent with 
each other. Moreover, if one needs to solve Laplace’s equation to find fringe 
field integrals, one might as well use the solution directly in a ray-tracing 
code. If one does go through this exercise, one discovers that the higher 
order aberrations are relatively insensitive to the ‘hardness’ of the quadrupole 
edges. This leads one to suspect that the aberrations are dominated by an 
intrinsic effect which has nothing to do with the detailed shape of the fringing 
field. Such is indeed the case. 

It turns out to be possible to find a canonical transformation which elim¬ 
inates the derivatives of V(z). I 11 our case, we wish to retain terms to 

4 th order in the Hamiltonian (3 rd order on force), and the transformation 
{x,p x ,y,p y ) —* (X, Px,Y, Py) has generating function 

G(x, P x , y, Py) = xP x + y P Y + ^(x 4 - y 4 ) + 

-Y( X *P x -y 3 P Y ). (6) 

0 

To the same order, this yields the transformation 

V , 

x = X + —X 3 
6 

Ho V' o 

P* = Px--X 2 P x + —X\ (7) 

The ^-transformation is obtained by replacing x,p x ,X, P x with y, p y , Y, P Y 
and V with —V. Note that outside the quadrupole, the transformed coordi¬ 
nates are the same as the original ones. 
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This yields the transformed Hamiltonian H*: 

H‘ = j(X 2 - Y 2 ) + \(P 2 X + P 2 ) + 

+ \(Px + Pyf - j(X 2 + r 2 )(Px - Pi) 

7V 2 V 2 

+ —(x 4 + y 4 ) - — x 2 y 2 . (8) 

24 v ’ 4 v ’ 

We can identify the terms: the first two are the usual linear ones; the third 
term is not related to the electric held (it is small and due to the fact that 
x' 7 ^ p x or, equivalently, tan 6 ^ sin 0); the 4 th term is also small and arises 
because a particle going through the quadrupole at an angle is inside the 
quad for slightly longer than one which remains on axis. See ref. ||3j for more 
complete physical derivation of the individual terms. 

The dominating higher order terms are the last two terms in eqn.[ 8 j Since 
there are no derivatives of V, we can directly write down the aberrations in 
the thin-lens limit: 

Ap ‘ = jk (h 3 ~ y' 2 ) ■ (9) 

with a similar expression for A p y . L and / are the quadrupole’s effective 
length and focal length. The fractional focal error is found by dividing by 
the linear part A 0 p x = —x/f: 


A f_x 
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for x, and similarly for y. 


2.2 Magnetic 

In magnetic fields, the canonical momentum p contains the vector potential 
A so that the time-based Hamiltonian is 


Hx = 


2 m 


p — qA 


( 11 ) 


As before, we use the invariant p 0 = \/2mH T to normalize the momenta, 
convert to z as independent variable, and expand the square root, keeping 
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terms up to 4 th order: 

1 r i 

H ~ — A z + — ( p x — A x ) 2 + ( p y — Ay ) 2 + 

+ g (Px — A x ) 2 + ( Py - Ay ) 2 . (12) 

To this order, the vector potential for quadrupole strength k(z) is 

A x = -jxy , Ay = —x 2 y (13) 

K = - t ^ 2 - y 2 ) + - y 4 )i 

and the Hamiltonian can be written: 

H = \ k (x 2 - y 2 ) - - y 4 ) + pI + pI + 

+ ypx - xp y ) + i (p 2 x + p 2 y f ■ (14) 

The generating function which will eliminate derivatives of k is 
G(x, P x , y, P Y ) = xP x + yP Y + ^(x 4 - y 4 ) + 

[(^ 3 + 3 xy 2 )P x - (3 x 2 y + y 3 )Py] , (15) 

which, to the same order yields transformation 

x = X + ^(X 3 + 3 XY 2 ) (16) 

p x = P x -^[(X 2 + Y 2 )P x - 2XYP y ] + ^X 3 , 
and similarly for (■ y,p y ). The transformed Hamiltonian is 
H- = |(A ' 2 - Y 2 ) + i(ff + P 2 ) + 

+ gf/’l + p ?) 2 - |(v 2 + r 2 )(Pl - Pp 

+ 1 )(A ' 4 + v 4 ) + yX 2 V 2 . ( 17 ) 
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Notice the similarity to eqn.[SJ in fact all terms are identical except the last 
two, which only differ in their coefficients. Applying the same procedure 
as in the electrostatic case, we write down the fractional change in focusing 
strength: 


Mx_ 
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3 Discussion 

Formulas [TU] and [TH] are handy for quickly evaluating the importance of 3 rd 
order aberration. They also show that for fixed focal length, the only way 
of reducing the aberration is by lengthening the quadrupole; the fraction of 
aperture used is not important; for a given effective length, the absolute size 
of the aperture is not important; the shape of the ends of the electrode is 
not important. 

Comparing the two formulas, we see that for roundish beams (x ~ y ), 
electrostatic and magnetic quads yield similar aberrations: they are in the 
ratio of ~ For cases where one transverse dimension is large compared 
with the other, and it is important to maintain the quality in the larger 
dimension, magnetic quads are better by a factor of However, for the 
more common case where it is more important to maintain the quality of the 
higher quality dimension, electrostatic quads win by a factor of 2. 

Results from using the above Hamiltonians are in agreement with those 
from using the commonly used codes GIOS and COSY, provided fringe field 
cards are used. In both of those codes it is possible to perform a 3 rd order 
calculation with quads which have no fringe fields. This gives incorrect and 
actually completely unphysical results. In essence, omitting the fringe field 
cards in those codes describes a situation where the particle traverses non- 
Maxwellian fields. For example, GIOS, since it does not use the scalar value 
of the potential field, does not obey conservation of energy when fringe field 
cards are omitted. 

The hard-edge case is correctly described in GIOS by including fringe field 
cards and setting the quadrupole aperture to zero, or, equivalently, setting 
all the fringe field integrals to zero. This is a useful approximation since 
the results are usefully close to reality and yet one needs not worry about 
specifying realistic fringe field integrals. This does not work in COSY, since a 
zero aperture forces an infinitesimal integration step-size. A better solution 
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would be to build in the hard-edge kicks and use these as default when no 
fringe held is specified. 

The required hard-edge kicks at the entrance to the quadrupole are de¬ 
rived directly from equations [7] and [16] The reason is that we know that the 
transformed coordinates (X, Px, Y ., Py ) do not experience any singular forces 
in the hard-edge limit. Therefore, the kicks for those coordinates are all zero. 
So the kicks for the untransformed (x,p x ,y,p y ) for the electrostatic case are, 


V 3 

= — X 


- 1 / 


-X Px 


-V 


1 / 


y 2 Py, 


Ax 

A p x 

Ay 

A Py 

and for the magnetic case are, 

Ax 
Ap x 
Ay 

Ap y = ^ [(x 2 + y 2 )p y - 2 xyp. 


^ 0 3 + 3 xy 2 ) 
k \x 2 + y 2 )p x - 2 xyp y 


-k 

~12 

k r / 


(3 x 2 y + y 3 ) 


(19) 


( 20 ) 


The kicks at the exit are, of course, opposite in sign. These agree with the 
GIOS case of zero fringe held integrals. See ref. [T]. 
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